Effect of anisotropy on universal transport in unconventional superconductors 
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We investigate the universal electronic transport for a mixed d x 2_ y 2+s-wa.ve superconductor in 
the presence of an anisotropic elliptical Fermi surface. Similar to the universal low-temperature 
transport predicted in a d 2 ,2_ y 2-wave superconductor with a circular Fermi surface, anisotropic uni- 
versal features are found in the low-temperature microwave conductivity, and thermal conductivity 
in the anisotropic system. The effects of anisotropy on the penetration depth, impurity induced T c 
suppression, and the zero-frequency density of states are also considered. While a small amount of 
anisotropy can lead to a strong suppression of the effective scattering rate and hence the density 
of states at zero frequency, experimental data suggests that large effects are restored by a negative 
s-component gap admixture. 



X 



PACS numbers: 78.30 .-j, 74.62.Dh, 74.25.Gz 
I. INTRODUCTION 

It is well known that in high-temperature supercon- 
ductors, owing to their layered nature and their reduced 
dimensionality, there exits significant anisotropy between 
various properties parallel and perpendicular to the cop- 
per oxide layers [|l|. In some materials such as YBCO, 
anisotropy is also seen between properties along the the 
a and b axes. These in-plane anisotropics are due to the 
presence of a CuO chain layer in addition to the CuC>2 
plane layers. Thus the overall band structure parallel to 
the layer is orthorhombic instead of tetragonal. On the 
other hand, the symmetry of the superconducting gap 
remains controversial in high-T c compounds. It is gener- 
ally believed to be d x 2_ y 2-wave, while many suggest that 
there is an additional s-wave component Hj. 

For a superconductor with order parameter with nodes 
on the Fermi surface, it is known that impurity scattering 
can lead to a finite density of quasiparticle states at zero 
energy (gapless excitations). As a consequence, while 
the effective impurity scattering rate 7 is quite different 
for different impurity concentrations and different scat- 
tering limits, the in-plane microwave conductivity satu- 
rates at low temperatures and is independent of 7 (or 
the impurity concentration). This universal transport 
phenomenon was first predicted by Lee ip| for a d x 2_ y 2- 
wave superconductor. Later Graf et al. pj found a similar 
universality in the thermal conductivity which is related 
to the microwave conductivity via the Wiedemann- Franz 
law. Recent thermal conductivity experiments done by 
Taillefer et al. || have confirmed these universal fea- 
tures. More recently, Wu and Carbotte || demonstrated 
how one can also study these universal features by do- 
ing channel-dependent Raman scattering experiments. It 
was found, in a d,^ ^2 -wave superconductor, that the 
low-temperature slope of the Raman intensity at zero 



frequency are universal in the A\ g and Big channels, but 
strongly dependent on the scattering rate (~ 7 2 ) in the 
Big channel |9|. 

The saturation of the low-temperature transport in an 
unconventional superconductor is a result of a cancella- 
tion between the value of the impurity-induced density 
of states at zero energy and the quasiparticle relaxation 
lifetime. It thus arises only if there exists nodes in the 
order parameter on the Fermi surface. In this paper, we 
consider the orthorhombic nature of the crystal structure 
in YBCO and explore whether an anisotropic universality 
exists in such a system. For simplicity, we study a system 
with an elliptical Fermi surface (assuming that the effec- 
tive masses are different along x and y directions). We 
have simultaneously considered a d x 2_ y 2-wawe gap with 
an additional s-wave component. We investigate vari- 
ous quantities, including the in-plane low-temperature 
microwave conductivity, thermal conductivity, and the 
London penetration depth. It will be shown that a small 
amount of anisotropy (the combined contribution from 
energy band anisotropy and gap symmetry anisotropy) 
can lead to a drastic change in the effective scattering 
rate in the superconducting state. As a consequence, the 
system exhibits an anisotropic universal feature in vari- 
ous transport quantities. More interestingly, experimen- 
tal data |l(]|ll|] seems suggesting that the effect of the 
large energy band anisotropy is largely compensated for 
by a negative s-wave component distortion to the d x 2_ y 2- 
wave gap. 

The content of this paper is as follows. In Sec. ||, we 
establish the Kubo formalism needed for subsequent cal- 
culations. In Sec. HI, we present our starting anisotropic 
model for a d x 2_ y 2+ s-wave superconductor with an or- 
thorhombic elliptical Fermi surface. In Sec. IV, the 
anisotropic universal microwave conductivity, anisotropic 
magnetic penetration depth, and anisotropic universal 
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thermal conductivity are presented. In Sec. [v[ we dis- 
cuss how the effective scattering rate and hence the den- 
sity of states at zero frequency are modified in such 
an anisotropic system. The effect of anisotropy on the 
impurity-induced T c suppression is also considered. A 



short conclusion is given in Sec. VI 



II. BASIC FORMALISM 

We derive first the formalism for the optical conduc- 
tivity in details. The microwave conductivity and the 
London penetration follow and are given by the real and 
imaginary parts of the optical conductivity in the limit 
of zero frequency. Later we will briefly derive the formal- 
ism for the thermal conductivity which shares a great 
similarity with that of the optical conductivity. 

The optical conductivity is given by Jl| 

<V,(0) = 
i 



O 



jr /ltf (q->0,ti/ B ->fi + ia)-M2(q-fO,0) , (1) 



where K^ v is the paramagnetic Kubo function given by 
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with Tr denoting a trace. In Eq. (|I|), the second term 

K$ is calculated in the normal state and with a minus 
sign corresponds to a diamagnetic response. The current 
vertex in Eq. (0) is given by 
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with Ti 's the Pauli matrices and £ k the electronic disper- 
sion relation of the superconducting layer. In Eq. (||), we 
have ignored the contribution to the vertex corrections 
due to the impurity potentials and superconducting two- 
particle pairing interactions. For isotropic impurity scat- 
tering, it is sufficient to use the bubble diagram at small 
q, while the inclusion of the pairing interaction vertex 
correction is negligible at the low frequencies of interest. 
However, the effect of impurities is fully included in the 
single-particle matrix Green's function G in Eq. (^). 

In terms of the particle-hole space, the single-particle 
matrix Green's function is given by 



G(k,iuj n ) 



iu) n Tp + gkT 3 + A k fi 

(iu) n ) 2 ~ I; 



(4) 



where uj n , £ k , and Ak are the impurity-renormalized 
Matsubara frequencies, electron energy spectrum, and 
gap. G is related to the noninteracting Green's func- 
tion Gq 1 (k,iuj n ) = iLu n To — £kT3 — A k Ti via the Dyson's 
equation 



G 1 (k,iw n ) = G 1 (k,iu) n ) - t(k,iui n ). 



(5) 



We shall solve the self-energy E due to the impurity 
scattering. By expanding E(k, ioJ n ) = E a (k, iuj n )f a 

at 

(a = 0, 1,3), one finds iu> n = iw n - E , Ck = Ck + E 3 , 
and A k = A k + Si. Employing the usual T-matrix ap- 
proximation, the self-energy is then given by £(k, iu> n ) — 
njT(k, k, iu) n ), where rii is the impurity density and 



T(k, k', iu> n ) =Vi(k,k')t 3 

+ ^(k, k")r 3 G(k", tw„)T(k", k' 



(6) 



Here Vi(k,k') = (k'\vi\k) is the impurity potential. If we 
consider only isotropic impurity scattering [i>i(k, k') = 
Vi], the T-matrix in (^|) is left only with frequency de- 
pendence and can be solved to get 



T(iuj n ) = 1 - ViT 3 G(iu> n ) 



ViT 3 



(7) 



with the integrated Green's function G(iui n ) = 
G(k, iu> n ). One can expand G(iu) n ) = G a (iu) n )f a 

a 

{a = 0,1,3) with G a {iuj n ) ee l/2£ k Tr[f Q G(k,iw„)]. 
For a superconductor with particle-hole symmetry, 
G^iujn) = and one obtains 



E (mj„) = 

Ei(iw n ) = 
E 3 (iw ri ) 



riiGo(iuJn) 



c 2 - G&iu n ) + G?(iw n ) 

-njGi(iuj n ) 

c 2 - Gg(iw n ) + G\{iw n ) 

CUi 

C 2 -G 2 ^„) + G 2 (^„)' 



(8) 



where c 2 = 1/vf. The effect of E3 is usually absorbed 
into the chemical potential and consequently £ k = £k or 
E 3 = 0. It is noted that for a system with a cylindrical 
Fermi surface and a pure d-wave gap, Gi = and hence 
Ei = which gives no renormalization effect on the gap 
(Ak = Ak). However, our present interest is to include 
this renormalization effect for an elliptical Fermi surface 
system and gap with a possible s-component appropriate 
to an orthorhombic system. The Green's function (|]) and 
above self-energies are solved self-consistently along with 
the gap equation 



Ak = ? E 5(k,k')Tr[f 1 G(k',zc„)], 



where g(k, k') is the pairing interaction. 
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III. ANISOTROPIC MODEL 

We consider a layered superconductor of a d x 2_ y 2+s- 
wave order parameter 
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A (k 2 x -kl + s), 



(10) 



where the constant s denotes the s-wave component 
which is small but can be positive or negative in general. 
For YBCO which is orthorhombic, the band structure is 
chosen to be elliptical 
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where ep is the Fermi energy and taking into account the 
CuO chain along the k y axis, the effective mass m x > m y . 
Our approach follows Kim and Nicol jl3| closely. It is 
convenient to use a transformation 



2m,, 



HI:, 



(A* = x,y) 



(12) 



and consequently Eqs. (|1C|) and (|11|) are transformed 
from the k frame to p frame according to 



Ao/(0) 



with 



cos(2</>) + a 
1 + a cos (2(f)) 



m x + m y 



(13) 



* > (14) 



and 

£(k x ,ky) -> £,(p x ,Py) 



m x + m y 



(^+p2)-^. (15) 



The denominator in (|14| ) was ignored in the work of Kim 
and Nicol jl3j. In (|13j), denotes the azimuthal angle 
in the p frame. The cylindrical Fermi surface case cor- 
responds to m x = my or a = 0. It is worth noting that 
the gap given in (O) viewed from the p frame is not 
a simple d+s- wave gap, A p = A (cos(2</>) + s), unless 
a = 0. Fig. [I] displays the gaps given by ([13]) and |L4| ) 
with various choice of a and s. The general feature of 
these pictures is that the gap nodes are shifted off the 
diagonals and for the case a — — s, the nodes are push 
back to the diagonal. These are to be compared by the 
pure d-wave case (a = s = 0) case. 

The momentum k sum can be transformed to p sum 
which in turn can be replaced by an integration 



/oo />2 
dt P / 
-oo JO 



2vr 



(16) 



with iV(0) = (m x + m y )/47rh defined as the density of 
states per spin on the Fermi surface in p frame. In the 
normal state in which A^ = and hence G\ — Ei = 0, 
the only nontrivial self-energy in (||) is E = n i Go/(c 2 — 
Gq). Using (|l~6|), one can easily work out the (isotropic) 
scattering rate «Eo = (l/2r), where in the Born scat- 
tering (c > 1) limit, l/2r = 2irN(0)nivf , while in the 
resonant scattering (c <C 1) limit, 1/2t = rii/2TrN(0). 
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FIG. 1. Order parameters on the Fermi surface given by 
Eq. (|l3|). (a): a = s — (pure d-wave gap); (b): a = 0.4, 
s = 0; (c): a = 0.4, s = 0.4; (d): a = 0.4, s = -0.4. The 
straight lines denote the nodal angles. 



Using the spectral representation for the imaginary fre- 
quency Green's function, analytically continuing to real 
frequency from imaginary frequency (iuj n — * to + iS), and 
then performing the frequency sum and momentum sum 
gives the real part of the Kubo function 



K »»(ty = — 7. — / ^7M<W^W 

2 Jo 27T 

x {[/(a;)-/(w-fi)]ReJ + _(n,w) 
- [f(u) + f(u-Q)]^Bl ++ (Q,u)} 

and the imaginary part of the Kubo function 

iV(0)e 2 



(17) 
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2 J Q 2tt 

x[f(u) - /(w - ft)] [Tm/ + _(n,w) +Im/ ++ (n,w)] . (18) 



In Eqs. (JTj) and (|18|), the vertices calculated on the 
Fermi surface are J x (<f>) — \f2epjrnx cos <j> and 7 y (</>) = 
we have defined 

Q'(Q + +Q') + A'(A + -A' + ) 



•v/2ej?/ my sin </> and we have defined 



.(n, w ) = -— 



.(ft,u;) 



(e+ + evK+ev 

d>L (£+ + + A'_(A + - A'_ 



(19) 



withcD± = iu> n (uj±id) ; tD± = iu> n (u>— ft±i#), A± = A(cj± 



i(5); A' ± = A(w - ft ± and £± 



sgn^^/u^. - Ai; 



£± = sgn(w — ft)y (<^±) 2 ~ (^±) 2 which are chosen 
to have branch cuts such that Im£ + , Im£^_ > and 
Im£_,Imfl < 0. The results of Eq. were ob- 

tained before by Hirschfeld et al. |l4j and many others. 
The microwave conductivity is given directly by K'' as 
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(20) 



f2->0 £1 

while the London penetration depth is given via K' as 

l 4?r r. 



with c the speed of light 



K' (Q = 0,T)-K'W(n = 0,T) 



, (21) 



IV. ANISOTROPIC UNIVERSAL TRANSPORTS 

A. Microwave Conductivity 

We consider the zero-temperature limit of the mi- 
crowave conductivity. When T — > 0, [f(oj) — f(u — 
f2)]/fi w df(oj)/dui w — as 17 — > 0. Consequently 
from (111) and (pPf), we have 



o-^(O) = iV(0)e 2 



JjMYjy (0)7 2 
[ 7 2 + A 2 (0)]i 



(22) 



where (• • •) denotes an average over the Fermi surface. 
In (||), 7 = -icj(w = 0) = i£ (^ = 0) and A(0) = 
A(cf),Lo = 0) = A o [/(0) + <5] with (5 = Si(w = 0)/A , 
respectively are the impurity-renormalizcd effective scat- 
tering rate and gap at zero frequency. For i « 1 and 
a + s < 1 of interest [see Eq. p4[) 1 , it is guaranteed that 
there exists nodes for the renormalized gap. The appear- 
ance of universal transport in unconventional supercon- 
ductors is intimately connected to the presence of nodes 
in the renormalized gap. When the anisotropics vanish 
(a = s = 0), f((f>) = cos(20) and (5 = which are cor- 
rect whatever concentration of impurity and correspond 
to the unrenormalized pure d-wave gap in a cylindrical 
Fermi surface case. We will discuss in the next section, 
the self-consistent results for 7 and 6 in the two different 
scattering limits, Born and resonant. 

In the case 7 = 7/A0 <C 1, the major contribution 
to the average in Eq. ( p2| ) comes from the small angular 
area around the nodes of the renormalized gap. One can 
then carry out the integration in (E2I) to get 



o-„„(0) 



4iV(0)e 2 
ttA 



iMohAM, (23) 



where 4>$ corresponds to the angle of nodes, i.e., f(<fio) = 
—8. The cross term cr xy (0) = because of the Fermi- 
surface average. More explicitly, Eq. (p3|) gives 



0-^(0) 



'yy 



(0) = 



2e 2 e F 7V(0) (1 + ar?) 2 (l + rj) 
TrA m x y/\ - ry 2 (l - a 2 ) ' 

2e 2 e F N(0) (1 + ar/) 2 (l - rj) 
irA m y ^/T~/f {I - a 2 )' 



(24) 



where we have defined r\ = cos(20o) = — {a + s + 5)/[l + 
a(s + 5)]. One can easily verify that when a = s = and 



hence 8 = corresponding to a d x 
cylindrical Fermi surface case (m x 

<Jxx{0) = Oyy{®) = 



_ y 2-wave gap in the 
m y ), 

(25) 

with n = kl/2-K for the electronic density. The universal 
feature of ( P5|) was first predicted by Lee raj. 

While the result in Eq. (§4) does not depend on 7, it 
is explicitly dependent of the gap shift 6 which in turn, 
is associated with the impurity concentration. In the 
large resulting anisotropy case [(a + s) 2 3> 1/(tA )] of 
primary interest (provided the impurity concentration is 
low), one can ignore 6 (compared to a and s) in Eq. (|24|), 
valid for both Born and unitary limits (see Sec. |v|). As 
a consequence, 



<7xx(0) 



'no 



(0) 



2e 2 e F iV(0) A x 
7rAo m 
2e 2 e F N(0) A y 



where we have defined 

1 — a 2 (1 — a — s 
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(1 + as) 2 VT 



or 



(27) 



p jyy ) is the plasma frequency in the x(y) 
6[. In comparison with experiment, the 
plasma frequency is a measured quantity and so <Ju(0) 
is a measure of An/ Aq. Note the extra factor of An in 
( p6| ) as compared with the well-known result (^5|) which 
is a(0) ^(fi 2 /4^)(lAA ). 



Eq. (|26j) represents the anisotropic universal feature 
in the microwave conductivity for an d^^^+s-wave su- 
perconductor with an orthorhombic band structure. The 
ratio 



cr ra (0) 



Q 2 A 

p,XX -T^XX 

O 2 ~A 

"p,j/j/ -"yy 



(28) 



a yy {0) 

should have some experimental consequence. We note 
that the low-temperature correction terms to an(0) will 
be proportional to T 2 /-f 2 with 7 the effective impurity 
scattering rate. Thus when T <C 7, the anisotropic uni- 
versal values can be attained. 

In the small anisotropy case [(a + s) 2 <C 1/(tAq)], 
while a and \s\ may be large individually, it can be shown 
that S ~ (a + s) ~ r\ are all small (see Sec. |y|). Conse- 
quently, one ignores r\ in (B4I) to obtain 



O"xx(0) 



,(0) 



n 2 



4^ 2 A (l-a 2 )' 

n 2 



4tt 2 A (1 - a 2 



(29) 



A ratio a xx {0) / a yy {0) 
limit. In the case of a 
and they reduce to Eq 



Q 2 

p,xx 



/n 2 

1 p,vv 



0, c7 X 40) 
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= (Xyy (0) in m 
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B. London Penetration Depth 

We consider here, the anisotropic London penetration 
depths within our anisotropic Fermi surface model to- 
gether with a d x 2_j / 2+s-wave order parameter. The re- 
sults are to be compared to the recent work of Kim and 
Nicol who effectively have considered the anisotropic 
penetration depth with a — because they ignore the 
denominator in (14). In the pure limit, using Eqs. (j^) 
and © it can be shown that the square of the inverse 
penetration depth is given by 



4- to = 1 

A,,,, 



\ " „,2 



k 



df(E k ) d/(e k ) 



(30) 



and is associated with the superfluid density. In the low- 
temperature limit, one obtains relatively simple results 
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fl 2 



n 2 

p.yy 



1 - 2 In 2 



2 In 2 



k B T 
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k B T 
Ao 



A, 



At, 



(31) 



where A xx and A yy are given in Eq. (p7|). 

Our result confirms Kim and Nicol's result when we set 
a = 0. The linear in temperature behavior in Eq. ( |3l| ) 
is expected for an unconventional superconductor with 
gap nodes on Fermi surface in the clean limit. When 
impurities are present, one instead expects a T 2 law for 
the low-temperature penetration depth (see Ref. p3|). 
In view of (|3l]), the discrepancy exhibited in a and 6-axis 
penetration depths can be studied in two ways. First, the 
zero-temperature anisotropic penetration depth ratio 



yy 



(0) 



n 2 

p. XX 



1 — a 



al(o) m„„ n, r i 

Secondly, the low-temperature slope ratio 
d[Xl x {0)/Xl x (T)]/dT\ 



IT^O 



d[XU0)/XUT)]/dT 



A,i 

X 



(32) 



(33) 
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To fit the experimental data of low-temperature pene- 
tration depth by Bonn et al. |l7j and the normal-state 
resistivity by Zhang et al. ]l8| ], it requires a — 0.4 in (|32] ) 
and then a value of s = —0.25 is obtained using ( |33| ) [p.9| . 

In Fig. ||, we plot the relative size of gap A(0) as a 
function of angle using parameters fit to low-temperature 
penetration depth slopes (alpha = 0.4, s = -0.25). As 
shown in Fig. the gap node is shifted from 45° and 
the overall gap magnitude is highly anisotropic. These 
may have observable consequence in future ARPES ex- 
periments. 
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FIG. 2. Relative magnitude of order parameter [given by 
Eq. (^)] vs. angle using parameters fit to penetration depth 
low-temperature slopes (a = 0.4 and s = —0.25). 
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^ 7m 0^7^(0) 



0,u>) +lml ++ (fl 



ou 

= 0,«)]. (34) 



As we seek the anisotropic universality in thermal con- 
ductivity, we are specifically interested in the T — > 
regime in which the inelastic scattering is unimportant. 
When T — > 0, the lo integration in Eq. (|34j) is limited to 
lu ss and consequently 



KjtuP -> 0) 



du) 



dx 



(35) 



Therefore the Wiedemann-Franz law (or Sommerfeld 
law) is satisfied. The above result is valid for both a and 
b directions. Previously Graf et al. have shown that 
for an isotropic Fermi surface band structure, the ther- 
mal conductivity saturates and obeys the Wiedemann- 
Franz law at low temperatures (T <C 7) for disordered 
apless superconductors. As emphasized by Graf et al. 
similar to the condition of universal microwave con- 
ductivity, the validity of the Wiedemann-Franz law at 
T —>■ is a unique feature of gapless superconductors 
in which the impurity induces a finite density of states 
at the zero energy. Our predicted anisotropic universal 
thermal conductivity can be written explicitly as [when 
( a + S ) 2 >l/(rA )] 



C. Thermal Conductivity 

In analogy to Eq. (|l^) which gives the microwave con- 
ductivity, the thermal conductivity is given by 
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Thus, similar to Eq. (| 

Kxx (0) 



,(0) 



^p,xx A xx 
TP A ' 

p,yy yy 



(37) 



Recently, an experiment done by Taillefer et al. ||] 
which measures the in-plane low-temperature thermal 
conductivity of YF^CusOg.g at different Zn substitu- 
tions for Cu has confirmed the universal feature in ther- 
mal conductivity. More recently, Chiao et al. [[ll| have 
observed the anisotropic universal thermal conductivity 
in YBCO. They found that the thermal conductivity sat- 
urates at low temperatures both along a and b directions 
from which an anisotropy k&(T — * 0)/n a (T — > 0) = 1.3 
is obtained filf . If we take a = 0.4 as predicted by 
anisotropic in-plane penetration depth ]l7fl and normal- 
state resistivity Jl^j data, it would imply a negative but 
large s = —0.60 in use of Eq. (|37|). Unfortunately such 
a large value of s reverses the order of the slopes in the 
penetration depth making the x direction steeper than 
the y direction in contradiction to experiment. For very 
small 1/t, it can also imply a very small value of 7 so 
that saturation effects will be seen only at very low tem- 
perature again in contradiction to experiment. On the 
other hand, with a = 0.4 and s = —0.60, the value of 
A xx ~ 1.9 makes it unnecessary to invoke a factor of 2 in 
the slope of the gap at the nodes needed in the analysis 
of Taillefer et al. Q who concluded that the actual gap 
grows out of zero a factor of 2 more slowly than in a pure 
d-wave case. 

The large negative value of s needed to explain the 
observed small anisotropy between a and b direction in 
the thermal conductivity is not believed to be physical. 
A more likely explanation may be in the observation that 
even for pure crystals of optimally doped YBCO, there is 
a large residual resistivity on the chains and none on the 
planes. This highly anisotropic residual resistivity has 
not been accounted for in our work. 



V. SELF-ENERGIES AND DENSITY OF STATES 

We discuss here the self-consistent results for the di- 
mensionless self-energies 7 and 8. One recalls that 7 = 
«Eo(w = 0)/Ao and 8 = £i(w = 0)/Aq. The self-energies 
Si in (||) are comprised of the integrated Green's func- 
tions Gi which are given by 



G (w) = -2i7rY(0) 



uj 2 — A 2 (0, uj) , 



and 



Gi(oj) = -2iirN(0) 



>Cu 2 - A 2 ((f>,uj) 



after the energy integration is done. For the case 8, 7 <C 1 
and |a + s| < 1 of interest (in which the gap is guaranteed 
to have nodes), we have found to leading order, 



G (w = 0) ~4iAT(0)7ln7 
Gi{uj = 0) ~ -AN(0)(a + s + S). 



A. Born Limit 



(40) 



Using Eq. 0) in (||), in the Born limit (c > 1), we 
obtain 



7 ~ e 



-7TT Aq 



ittAq 



(41) 



Typically tAq ~ 100 (providing the impurity concentra- 
tion is low), one thus retains 7 <§; 1 and 8 <C \a + s\ < 1. 
The latter validates Eqs. <M) and (ph. 



B. Unitary Limit 

In the unitary limit (c -C 1) of primary interest here, 
and when the resulting anisotropy is large such that (a + 
s) 2 > 1/(tA ), we found 



-4rA (a+s) 



4rA (a + s) 



(42) 



Clearly Eq. ( [i"2"| ) gives 7, S <C 1. Thus one can ignore S in 
@ to get Eqs. @ and @. Eq. @ indicates that the 
effective scattering rate is strongly suppressed compared 
to isotropic case. When the resulting anisotropy is small 
[(a + s) 2 < l/(rA )], though a and \s\ may be large 
individually, the system undergoes a crossover from an 
anisotropic regime to a "quasi-isotropic" regime in which 
the gap shift 8 is small and is linear in and about the size 
of (a + s). When a — —s, 8 vanishes as for a pure d- 
wave superconductor in an isotropic Fermi surface. More 
explicitly, to second order in (a + s), we find 



7 
8: 



At Ao In 70 

a + s 
1 + hi7o ' 



2rA ln7o(a + .s) 2 



7r(ln7o + l) 5 



(43) 



where 70 satisfies 7g In 70 



-it / (At A ). Typically 



(38) hi7o < —2 and with (a + s) 2 <^ l/(rA ), the second 
term in the expression 7 in Eq. ( ^3| ) is small. One thus 
obtains a 7 value similar to that for the case of a pure 
d-wetve superconductor with an isotropic Fermi surface. 

One should keep in mind that while when (a + s) 2 <C 
1/(tAo), the results for 7 and 8 are similar to those of an 

(39) 

isotropic system, the universal microwave conductivities 
are still anisotropic, as given by Eq. (f29|). They reduce 
to the universal value of Eq. (E5h only when a = 0. 
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FIG. 3. Plot of zero-frequency density of states N(0, a) 
scaled to N(0), in the unitary limit, as a function of a for 
different choices of scattering rate 1/AoT. Here we choose 
s = 0. 



C. Density of States 

Using the result of (|40|), the impurity renormalized 
density of states at zero frequency is given by 



N(0) 



2tt 



Go(w = 0) 



2N(0) 



7 In 7, 



(44) 



which is dependent on the impurity concentration (or 
1/t) and the magnitudes of a and s. In Fig. ||, we 
compute and plot the a dependence of N(0), in the 
unitary scattering limit, for three different values of 
1/(tAo), taking s = 0. Similar plot is made in Fig. |] 
using s = —0.25. As already predicted for the effec- 
tive scattering rate 7, the density of states is strongly 
suppressed when (a + s) 2 > 1/(tAq). For example, 
in case of 1/(tA ) ~ 0.01, N(0) is suppressed by 50% 
when (a + s) > 0.1 compared to (a + s) = case (see 
Figs. |3| and In comparison with experimental re- 
sults, the large discrepancy between a and 6-axis pen- 
etration depths [0 and normal-state resistivity |Q sug- 
gest that a large a ~ 0.4 is required so that N(0) is 
strongly suppressed. On the other hand, the low but 
finite-temperature universal observed value of the ther- 
mal conductivity || suggests that the effective scattering 
rate or density of states at zero frequency is not strongly 
suppressed. For one to be consistent with the other, it 
implies that a large a is accompanied by a large but neg- 
ative s which reduces the effect of a and leads to less 
suppression of the density of states at zero frequency. 



D. T c Suppression 

In the standard approach, we choose the pairing inter- 
action to have a separable form <?(k, k') = ~gf{<P)f((f>') 
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FIG. 4. Similar plot as Fig. | with s = -0.25. 

(g > 0), where f((f>) is given in (]l4|). This reduces the 
gap equation of Eq. @ to the form 



- = 27riV(0)T^ 

9 ,., 



V, /■(</)) A(0,i^)/A o 



+ A 2 (</>,w„), 



(45) 



where the prime denotes a cutoff. When T — > T c , it is 
sufficient to use the linearized version of (|4^), i.e., ignor- 
ing the A term in denominator. In a parallel manner, 
one can neglect the G\ term in the denominators of the 
self-energies in (|J). This enables one to obtain, valid for 
both Born and unitary limits, 



A(0,w„) = A 



1 



1 </(</>)> 



2r 



2t \u n 
sgn(w n ). 



(46) 



The second line of ( |46|) is identical to the solution for the 
normal state, as expected. Thus using (Q), the linearized 
version of mH) is reduced to 



- = 4nN(0)T V 
q 



1 



2tui, 



-</(£» 2 



n>0 



1 

27 



(47) 



We have found quite generally 

t c \ </ 2 w>-</w> 2 



In 



T, 



r() 



(PW) 



l/2r 



(48) 



where T c q is critical temperature in the pure case. 
In our case [see (pi)], 

i-vT 



(/ 2 W> 



Of 



(1 + 2as) + s 2 



1 - VI - a 2 



(49) 



7 



and to leading order, ((/ 2 ) - (/> 2 )/(/ 2 ) = (2 - a 2 )/[2 + 
As (a + s)]. When l/2r <§; T c q, we have obtained using 



Tc_ 



1 



2 + 4s(a + s) 



8rT f 



(50) 



which depends linearly on the impurity concentration. 
When l/2r > T c0 , we obtain 



T 

± c 

T c o 



1 



4:TttT c0 



(51) 



It is shown clearly in Eqs. (|5(j) and ( [5l| ) that, in contrast 
to the cylindrical Fermi surface and pure d-wave gap case 
(a = s = 0) in which T c is strongly suppressed compared 
to T c o, the presence of anisotropy (finite a or s case) 
reduces the effect of the suppression of T c . 



VI. CONCLUSIONS 

In this paper, we have studied the low-temperature 
microwave conductivity and thermal conductivity for a 
d x 2_ y 2+s-wa,ve superconductor with an orthorhombic el- 
liptical Fermi surface. Similar to the universal behav- 
iors found in the microwave conductivity Q and thermal 
conductivity for a d x 2_ y 2-w&ve superconductor with a 
cylindrical Fermi surface, anisotropic universal features 
are found in the present case. The effects of Fermi sur- 
face orthorhombicity and additional s component to the 
gap on the penetration depth, impurity induced T c sup- 
pression, and the zero- frequency density of states are also 
considered. It is found that, compared to the cylindri- 
cal Fermi surface and pure d x 2_ y 2-wewe gap case, a small 
amount of anisotropy (either band anisotropy or gap ad- 
mixture) will lead to a strong suppression of effective 
scattering rate and thus the density of states at zero fre- 
quency. Nevertheless, experimental data suggests that a 
large band structure anisotropy effect is compensated by 
a large but negative s-wave gap component. 
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